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Abstract 

This paper discusses the stabilization system 
analysis and orbital performance of the GEOS-A 
gravity-gradient satellite which was launched in- 
to orbit on 6 November 19<?5* The analysis is 
performed with the aid of a rather general digital 
simulation of the non-linear differential equa- 
tions of motion. It includes the effects of 
gravity -gradient forces, orbital eccentricity, 
eddy current dampers, residual magnetic dipoles, 
solar radiation pressure, and thermal distortion 
of the spacecraft's stabilization boom. 

The attitude measuring system aboard the space- 
craft is described, and processing of the attitude 
data is discussed. Description of failures in the 
attitude determination system and discussion of 
the subsequent effects on the experimental atti- 
tude results are included. Experimental attitude 
results are compared with the results of the digi- 
tal simulation. General agreement is obtained. 


Nomenclature 

a = s emi -ma j or axis 

E = light intensity 

H - magnitude of geomagnetic field 

H q = magnitude of geomagnetic field at mag- 

netic equator 

I , I ,1 = inertia of spacecraft about body fixed 

^ z principal axes 

X ,1 ,1 = inertia of spacecraft about pitch, roll 

^ r ^ and yaw axes 

k = gravitation constant 

M = mean anomaly 


n 

r 

t 

€ 


= mean motion 

= radial distance from cm of earth to 
cm spacecraft 
= time 

- orbit eccentricity 


Introduction 

The purpose of this paper is to discuss the de- 
sign analysis and flight performance of the stabi- 
lization system of one of the Geodetic Earth Orbit- 
ing Satellites, GEOS-A. Geocentric stabilization 
is achieved in a completely passive manner by 
utilization of the gravity -gradient principle and 
a passive energy dissipator. The system analysis 
is accomplished through a digital solution- of the 
equations of motion which describes the orienta- 
tion of the spacecraft in orbit. Flight perform- 
ance is obtained by an on-board attitude measuring 
system which determines the orientation of the 
satellite -sun line and magnetic field vector in 
the satellite reference system. A comparison of 
the theoretical and measured orientation of the 
spacecraft concludes the paper. 


The repeated success of passive gravity -gradient 
stabilization systems at low altitudes has foster- 
ed interest in the use of this technique for 
communication and meteorological satellites at 
higher altitudes with correspondingly more severe 
requirements on system performance. Since the 
gravity-gradient restoring torque is inversely 
proportional to the cube of the distance of the 
spacecraft from the center of the earth, increased 
altitude requirements demand a more precise knowl- 
edge of the effects of the disturbing forces. At 
synchronous altitudes the restoring torque is two 
orders of magnitude smaller than at low altitudes. 
For this reason, verification of the performance 
of gravity- gradient stabilization systems becomes 
an important consideration. 

GEOS-A was conceived as part of the National 
Geodetic Satellite Program to study geometric and 
gravimetric geodesy. Geometric geodesy utilizes 
the principles of triangulation to map the surface 
of the earth. The study of the gravitational field 
of the earth is the objective of gravimetric geod- 
esy. Five geodetic measurement systems are employ- 
ed to provide simultaneous evaluation of their 
performance and to ensure highly accurate measure- 
ments. The five geodetic systems are: 

1. A radio Doppler system consisting of three 
coherent frequencies for Doppler- shift determina- 
tion of the precise orbit. 

2. A radio ranging system which receives and 
retransmits ground based radio signals sent in 
sequence from four stations. The range is meas- 
ured by observing the total phase shift of an accu- 
rately known modulation frequency. 

3. A range and range rate system which is used 
to determine the changing range and radial velocity 
of the satellite by measuring the phase difference 
and doppler shift of a signal transmitted from a 
ground station and retransmitted by the satellite. 

4. Corner cube quartz laser reflectors to deter- 
mine the position of the satellite by measuring 
the time increment required for the laser beam to 
travel from the station to the satellite and back 
to the station. 

5- An optical beacon system to be photographed 
against a background of stars to provide angle 
measurements. 

GEOS-A is the largest and most' complete satellite 
ever constructed exclusively for the study of 
geodesy. It was designed and built for the Nation- 
al Aeronautics and Space Administration by The 
Johns Hopkins University Applied Physics Laboratory. 
Scientific and engineering direction of the project 
was provided by Dr. George C. Weiffenbach the pro- 
ject scientist and Ralph E. Willison the project 
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engineer. An artists conception of the spacecraft 
is presented in Fig. 1. It weighs 386.6 lbs. and 
is an eight -sided truncated pyramid. Each of the 
l6 flat surfaces carries a panel of solar cells to 
convert solar energy into electricity. 

For the GEOS -A mission an earth pointing satel- 
lite is desirable because it makes optimum use of 
the radiated and reflected optical and radio power. 
The most severe requirement of the stabilization 
system was imposed by the optical beam pattern. 

Hi is system was designed to compensate for atmos- 
pheric attenuation with large zenith angles on the 
assumption of achieving earth pointing accuracy to 
within 5 degrees. Fig. 2 illustrates the optical 
beam pattern of the lamp- reflector combination vs. 
the angle off the longitudinal axis of the space- 
craft. A comparison with the optical pattern of 
the ANNA 13 satellite is presented. 

Attitude Determination System 

The attitude determination problem is defined 
mathematically as the evaluation of a transforma- 
tion matrix which relates the orientation of the 
spacecraft to some specified reference coordinate 
system. The nine elements of this matrix are 
direction cosines and, of course, are dependent; 
they can be expressed in terms of three independ- 
ent variables. The transformation matrix can be 
determined as a function of time if two vectors 
can be specified simultaneously in both the satel- 
lite and reference coordinate systems. These vec- 
tor components, six scalar quantities in all, will 
provide redundant information for the evaluation 
of the three independent variables of the transfor- 
mation matrix so that in-flight calibration of the 
attitude system and statistical estimation of the 
attitude is possible. Two rather well defined 
vectors which can be measured in the spacecraft 
with a minimum expenditure of weight and power are 
the unit vector of the satellite -sun line and the 
vector geomagnetic field. On GEOS-A the measure- 
ment of these two vectors forms the basis of the 
attitude detection system. 

The instrument to measure the geomagnetic field 
components consists of three second-harmonic type 
flux-gate magnetometers. Calibration was perform- 
ed at the Naval Ordnance Laboratory 1 . In addition, 
partial calibration of the instrument can be achiev- 
ed in-flight by superposition of known magnetic 
biases on the three detectors sampling the ambient 
magnetic field. 

The satellite-sun unit vector components are 
measured by an array of specially designed detec- 
tors utilizing analog solar cells. After consider- 
able study at APL a system consisting of three 
types of detectors evolved^»3. Each sensor module 
is designed to contain: l) an analog cosine sun 
sensor having a response grossly approximating a 
cosine curve as a function of the angle of inci- 
dence; 2) a detector which has a linear response in 
the vicinity of normal incidence; 3) a combination 
of two detectors providing a step response at each 
of two angles in the vicinity of normal incidence, 
see Fig. 3* Flight hardware, selected on the basis 
of similarity of response characteristics, is 
calibrated before mounting on a spacecraft. 


gravity-gradient boom of the spacecraft, hereafter 
referred to as the +Z direction, and the remaining 
four modules were mounted around the equator of 
the spacecraft pointing in mutually orthogonal 
directions denoted hereafter as + X, - X, + Y, and 
- Y. No module was mounted so as to point in the 
-Z direction since in the normal gravity- gradient 
stabilized orientation the - Z axis is directed 
earthward. To minimize the effects of reflected 
light from the earth, shades were used below the 
equatorially mounted modules. 


Experience has shown that the calibration con- 
stants for analog cosine detectors change during 
flight. Black has developed a program for recali- 
brating these instruments in f light J+ . After re- 
calibration , attitude data is processed by the 
Satellite Attitude Determination program (SAD), 
also developed by Black5>6. This program requires 
the orthogonal vector components of the earth* s 
magnetic field vector and the orthogonal vector 
components of the satellite-sun unit vector, both 
measured in the satellite* s coordinate system. 


Stabilization System Description 

Gravity-gradient stabilization utilizes the 
interaction between the inertia ellipsoid of the 
spacecraft and the gravitational field of the earth 
to provide geocentric orientation. Initial mis- 
orientation and perturbing torques induce satellite 
motion about its equilibrium position so that a 
damping scheme must be provided to dissipate this 
energy. Passive stabilization and damping offer 
the significant advantages of high reliability 
and a long lifetime . 

The GEOS-A stabilization system consists of a 
single motorized extendible boom manufactured by 
the DeHavilland Company; at its end is a General 
Electric magnetically anchored eddy current 
damper^. The boom motor is controlled by ground 
command and backed up by an automatic stop switch 
activated by detents on the boom tape. This sys- 
tem provides bi- stable equilibrium about any axis 
normal to the local vertical. As desired for 
thermal uniformity there is no preferred orienta- 
tion about the longitudinal axis. 

A preliminary analysis indicated that this sys- 
tem, which would weigh approximately 11 lbs, 
would meet two stabilization requirements: l) oper- 
ational status within 15 days after launch, and 
2) a maximum steady-state libration amplitude of 
5 degrees. An important consideration in the 
selection of this system was that it had been 
successfully demonstrated in orbit. The emphasis 
on this factor was consistent with the mission, 
in that GEOS-A is basically a geodetic satellite 
not a gravity-gradient experiment. 

The launch vehicle procured was the new improved 
Delta in which the last stage is spin stabilized. 
When separation occurs, the spacecraft rotates at 
a rate of approximately 150 rpm. A yo-yo despin 
device reduces this rate to about 1 to 3 rpm. As 
a despin backup, three orthogonal ferromagnetic 
rods encased in aluminum jackets are provided-. 

These act as ' an eddy current despin device at high 
spin rates. 


Five solar detector modules were mounted on the 
GEOS-A spececraft. One module was mounted so as 
to be pointing in a direction parallel to the 


Proper stabilization of the satellite can be in- 
sured if the extension of the boom is initiated 
when the boom axis is pointing upwards and within 
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55 degrees of the local vertical. However, there 
was no assurance that this orientation would be 
obtained during a pass over the command station at 
the Applied Physics Laboratory. Another factor was 
that hazardous temperatures could be reached by 
the spacecraft which was thermally designed to 
operate as a gravity- gradient satellite. Both of 
these considerations indicated that boom extension 
might need to be initiated when the boom axis is 
in a non-preferred direction possibly resulting in 
upside-down capture. Inversion of 'the spacecraft 
could then be accomplished by a properly timed 
sequence of boom retraction and extension. 

The eddy current damper designed and built by 
the General Electric Company consists of a magnet 
assembly which is free to lock onto the magnetic 
field of the earth. Rotation of the magnet assem- 
bly relative to a conducting spherical shell 
induces eddy currents in the conductor which 
impedes the motion of the spacecraft. To prevent 
physical contact between the magnet and the con- 
ducting sphere a diamagnetic shell of pyrolytic 
graphite is used as a centering device. 

Attitude Simulation 

The analytical theory upon which passive gravity- 
gradient stabilization is based is well known. The 
analyses, however, tend to be complicated by the 
non-linearity and coupling of the equations of 
motion and are intractable without simplifying 
assumptions about the perturbing forces. The three 
non-linear differential equations for the attitude 
motion of a spacecraft in orbit about a gravitat- 
ing mass were developed in terms of Lagrangian 
mechanics®. The center of mass of the satellite is 
constrained to an orbit of modest eccentricity 
having a longitude of the ascending node and an 
argument of perigee which are linear functions of 
time. These secular variations result from the 
even zonal harmonics of the gravitational field, 
principally the Jr, term. The coupling between 
orbital and body motions of the spacecraft is neg- 
lected. The attitude is expressed in terms of a 
general set of Euler angles because equations of 
motion in terms of any one set of these variables 
always have a singularity. This generality allows 
the use of a dual set of angles to specify the 
attitude without restrictions. 

A rather general digital computer program titled 
Digital Attitude Simulator (DAS) has been developed 
which' is applicable to a wide class of satellite 
configurations. It includes the effe'cts of 
gravity- gradient torques, residual dipoles, solar 
radiation pressure, thermal distortion of gravity- 
gradient booms and passive dampers of the eddy 
current and hysteretic types. A high order expan- 
sion of the rotating geomagnetic field in terms of 
spherical harmonics is included. The integration 
of the non-linear second-order differential equa- 
tions of motion over the time domain is accomplish- 
ed by a fourth -order Runge Kutta integration. 

The eddy current damper is simulated on the 
assumption that its magnet assembly is locked to 
the geomagnetic field. This is valid for the rath- 
er slow rates of motion of the spacecraft relative 
to the magnetic field and the large magnetic mo- 
ments of the magnet assembly which are used. The 
magnitude of the damping torque is proportional to 
the angular velocity of the conducting shell rel- 
ative to the magnet assembly. The thermal distor- 


tion of the -gravity-gradient booms is handled in a 
quasi- static manner. Each boom is assumed to bend 
with a constant radius of curvature in the plane 
containing the undeformed boom and the sun. The 
boom deformation results in changes in the mass 
distribution parameters. The radiation pressure is 
computed for the booms in their undeformed position. 
That part of the incident radiation which is 
reflected is a function of the surface character- 
istics and is in general assumed to be of both the 
specular and diffuse types. The effects of aero- 
dynamic torques are negligible at the altitudes of 
interest. 

A modification of the general program is used to 
analyze inversion maneuvers by taking into account 
the time dependent mass properties of the space- 
craft. It is assumed that the boom retraction and 
extension rates are different and are polynomial 
functions of time. 

Pre launch Design Study 

Pre launch design studies were carried out inde- 
pendently by the Applied Physics Laboratory and 
the General Electric Company. In general, there 
was very good agreement. At APL the design study 
was performed by utilizing the general digital 
simulation described above. Por this particular 
spacecraft and orbit the major perturbing torques 
arise from the gravitational, magnetic, and solar 
environments. The gravitational field >through the 
orbit eccentricity, induces a steady-state pitch 
libration. A linear theory in Appendix A gives an 
approximation to the amplitude and frequency of the 
perturbation. The magnetic torques arise from the 
interaction of the damper and residual dipole with 
the geomagnetic field. The coefficients used in the 
spherical harmonic expansion of the field consti- 
tute the Jensen-Cain model which has an rms 
accuracy of 9 percent for order 2 (8-terms) and 1 
percent for order 4 (24- terms )9. 

Studies of the transient and steady- state perform- 
ance of the spacecraft were made for the nominal 
orbit which had a perigee of 600 nm and an apogee 
of 800 nm. The eccentric orbit, which is undesir- 
able for accurate stabilization, is necessary for 
accurate gravimetric geodesy for which the perigee 
need be well defined. The nominal orbit parameters 
and mass properties of the spacecraft are given in 
Table I. 

The requirement on the transient performance of 
the stabilization system was that the spacecraft be 
in operational status, that is correctly oriented 
and damped to near steady-state, within 15 days 
after launch. This 15- day interval includes despin 
from about 1 to 3 rpra to near orbital rate, waiting 
foif favorable orientation of the longitudinal axis, 
erection of the stabilizing boom, damping to steady- 
state and, if required, the inversion maneuver and 
subsequent damping. 

Approximations to the time constants of the de- 
spin and libration motions are possible. During the 
despin the large angular velocities of the space- 
craft Induce rotational motion of the magnet 
assembly so it is not locked to the magnetic field. 
The simple exponential law for eddy current damp- 
ing which gives a time constant proportional to 
i/c is optimistic. The amplitude of the libration 
motion after boom extension depends primarily on 
the initial orientation of the longitudinal axis. 
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The preferred orientation is that the boom axis is 
pointing up and near parallel to the local vertical. 
As the loom is extended, the work done by the gravi- 
tational forces will be small so there is a tenden- 
cy to conserve angular momentum. The large inertia 
change will reduce the rate of rotation about all 
axes except about the longitudinal to near zero. 
Since the equilibrium orientation of the gravity- 
gradient satellite has an in-plane angular velocity 
of one cycle per orbital revolution, pitch libra- 
tions will result. For this reason, a study of 
the pitch motion best illustrates the effectiveness 
of the damper. A linearized solution for the trans- 
ient pitch lib rat ion applicable to polar and near 
polar orbits is given in Appendix B. The time con- 
stant of the envelope of the transient response is 
proportional to I /c, as was the despin time con- 
stant. This solut¥on also shows that, as expected, 
the damper induces a disturbing torque which is in- 
versely proportional to Ip/c. Any improvement in 
the transient performance c>y decreasing Ip/c will 
result in a larger disturbing torque and poorer 
steady-state performance. Preliminary estimates 
indicated that for boom lengths of 50 to 60 feet 
damping coefficients of 60,000 to 80,000 dyne- era- 
se c would satisfy. the transient requirements. Pre- 
liminary simulation results indicated that this 
configuration would also satisfy the steady- state 
requirements. 

A contract was then let to the General Electric 
Company for the fabrication of the damper and 
analysis of the stabilization system. 

Concurrent with the fabrication of the satellite, 
a detailed simulation study of the steady- state 
performance was initiated. For GEOS-A the steady- 
state lib rat ions result principally from the orbit 
eccentricity, magnetic damper, residual magnetic 
dipole, radiation pressure and thermal distortion 
of the gravity- gradient boom. It was recognized 
that the magnet assembly of the damper, which was 
not caged during launch could induce a residual 
magnetic moment in the batteries which are adjac- 
ent to it. Experiments confirmed this and pro- 
vided bounds for the magnitude of the residual 
magnetic moment of the spacecraft^. This also 
meant that the magnetically anchored damper would 
be ineffective during the despin phase since its 
magnetic assembly would lock to the batteries. 

This problem could be circumvented by extending 
the boom a small distance. Once the characteristics 
of the damper had been finalized, by the General 
Electric Company, to a weight of 7 lbs* a damping 
coefficient of 70,000 dyne-cm-sec and a pole 
strength of 25,000 pole-cm, the remaining free 
variable in the steady- state analysis was the boom 
length. An increase in boom length, which increases 
the inertia of the system, tends to minimize the 
effects of the residual dipole and eddy current 
damper but increases the effects of radiation pres- 
sure and thermal bending. Therefore some optimum 
boom length exists which would minimize the total 
steady-state amplitude. Preliminary results indicat- 
ed that a boom length of 50 ft would meet the 5 
degree requirement. The boom tape was cut at 56 
feet to provide a margin of 6 feet for any contin- 
gency, and a detent made in the tape at the 50 ft 
length to provide automatic motor shut off. 

Subsequent analysis by the General Electric 
Company and APL indicated that the 50 ft length 
would not meet the five degree requirement over 


some portions of the spacecraft’s lifetime. The 
amplitude of the solar induced libration is a 
function of the relative positions of the orbit 
and the sun. A brief analysis of the twice orbit- 
al roll component of the solar forcing function 
indicated that the maximum roll amplitude is a 
function of the angle between the earth-sun line 
and line of nodes and the amount of the orbit 
that is shadowed. The analysis indicated that for 
the nominal orbit the maximum static roll ampli- 
tude induced by thermal -bending occurs when this 
angle is 65 degrees and when perigee is in the 
shade, see Fig. 4. Digital simulation results, 
also given in Fig. confirmed this analytical 
theory and indicated that for perigee in the shad- 
ow a maximum roll amplitude of 8.5 degrees was 
possible for a 50 ft boom. Simulation results for 
40, 50, and 60 ft boom lengths are given in Fig. 5 . 
These results indicated that the 40 ft boom would 
satisfy the steady-state requirements under the 
worst possible case of solar aspect and location 
of perigee . Since the maximum amplitude appeared 
not to be very sensitive to small errors in boom 
length, it was decided that the proper length 
could be obtained by timing the boom extension 
and counting the number of turns of a gear in the 
motor drive . 

Estimates of the coast time for the inversion 
maneuver indicated that the total time required 
to retract the bocm, coast, and extend the boom 
precluded the operation during a pass over a sin- 
gle station. However, the maneuver could be per- 
formed over an arc of the orbit covered by two 
ground stations whose coverage only slightly over- 
lapped. The NASA STADAN Station at the University 
of Alaska in Fairbanks together with the Applied 
Physics Laboratory station at Howard County, 
Maryland met this requirement . 

Post -Launch Analysis 

GEOS-A was successfully launched into orbit at 
1836 TJT on day 310 in 1965 from the Eastern Test 
Range . The spacecraft was officially designated 
1965 89 A and titled by the National Aeronautics 
and Space Administration — Explorer XXIX. While 
the boost phase of the launch was normal, the sec- 
ond stage apparently burned to depletion. This 
resulted in obtaining the target perigee of 600 
nautical miles but an apogee of 1230 nautical miles 
which is 4-30 nautical miles higher than planned. 

On the first pass over Alaska the spin rate was 
observed to be about one rpm. This confirmed the 
successful deployment of the yo-yo despin device . 
No reduction in spin rate was observed during the 
next 14 hours indicating that the magnetically 
anchored eddy current damper was indeed locked onto 
the batteries of the spacecraft. 

The extension of the boom to its 3 foot detent 
was commanded from the Applied Physics Laboratory 
during the beginning of the 1212 UT pass on Day 
311. This operation was completely successful. By 
observing the decay rate over several passes a 
time constant of 100 minutes was determined. By 
1900 ITT on the same day the spin rate was observed 
to be essentially zero. After despin, the spacecraft 
was observed, through telemetered magnetometer 
data,to maintain a near fixed orientation with re- 
spect to the geomagnetic field. The residual mag- 
netic dipole induced by the damper resulted in the 
spacecraft being near magnetically stabilized. Re- 
duction of the magnetometer data, by Mobley, 
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resulted in the direction cosines of the residual 
magnetic moment vector in the spacecraft’s frame 
of reference-^-. 

As described above, the plan to stabilize the 
spacecraft was to extend the bocm to 40 feet when 
either the boom axis achieved an orientation which 
would ensure proper capture or the battery tempera- 
tures approached dangerous levels . On the 805 Iff 
pass on day 312 over the Applied Physics Lab or at cry 
the battery temperatures were observed to be below 
the lower flight acceptance limit so the boom was 
extended. The motor was started at 819 Iff and stop- 
ped at 827 Iff for an elapsed time of 460 seconds . 
This time interval corresponded to the required 
time for the boom to extend to about 40 feet when 
tested on the ground. The boom length was also 
inferred from telemetered data indicating the 
number of revolutions of a gear in the motor drive. 
This method indicated that the boom length was 
36.8 ±1 ft. On subsequent passes, the temperature 
increased to nominal values and magnetometer data 
confirmed that gravity-gradient capture was achiev- 
ed. Unfortunately the spacecraft was inverted. 

Since the solar detectors had earth shades, 
these earth shades now became sun shades and pre- 
vented the solar detectors frcan being fully illu- 
minated. Consequently, no definitive experimental 
attitude data could be obtained. 

At this time a simulation study of the inversion 
maneuver was undertaken utilizing the best esti- 
mate of the actual boom length and tracked orbit. 
The purpose of the study was to determine the 
coast time interval which would maximize the prob- 
ability of success of the maneuver. Station alerts 
based on the tracked orbit indicated the times 
favorable for the two -station maneuver. A deci- 
sion was made to attempt it at 1642 Iff on day 3^7 • 
The 5 day interval frcan time of gravity-gradient 
capture allowed most of the transient motion to be 
damped out. The sequence of events in the maneuver 
was: to retract the boom to its 15 ft detent, and 
after the coast phase, to reextend the boom. It 
was decided to reextend the boom to the 50 ft de- 
tent in case a second inversion maneuver would be 
required. The opt imum boom length could also be 
obtained more accurately at sane future time by 
retraction from the 50 ft detent rather than meas- 
ure the length by timing and counting gear revolu- 
tions during the inversion maneuver. 

Different quadratic bocm extension and retrac- 
tion rates obtained frcm ground experiments with 
the GEOS -A boon were used in the simulation. The 
major uncertainty besides initial boom length was 
in the pitch librational motion at the beginning 
of the boom retraction phase. Since the major 
pitch disturbing torque arises frcan the orbit ec- 
centricity, a good approximation to the pitch mo- 
tion could be obtained frcan the linear theory 
given In Appendix A. The results of the simulation 
study are given in Fig. 6 . The nominal initial con- 
ditions referred to are those computed from Eq. 

(A. 8 ) at the mean anomaly corresponding to 1642 Iff 
on day 317 . Small errors in the initial pitch rate 
have a much greater effect on the optimum coast 
time than small errors in pitch orientation. To 
study this phenomena, runs were made in which 
pitch rates corresponding to pitch librations of 
±3 and ±5 degrees were added to the nominal initial 
conditions. A coast time of about 19-25 minutes 
was selected as that which would give the highest 


probability of success . 

The inversion maneuver was initiated at the NASA 
STADAN station in Alaska and completed over the AIL 
station. The boom retraction and extension phases 
took 295 and 48l seconds respectively which com- 
pare to the ground test values of 282 and 485 sec- 
onds which were used in the simulation. The differ- 
ence in boom retraction t ime indicated that the 
bocm length could have been at the upper limit of 
the estimated length. Although lack of experimen- 
tal attitude data, larger than nominal eccentric- 
ity, and uncertainty in boom length combined to 
minimize the probability of success of the inver- 
sion maneuver, surveillance of the attitude deter- 
mination system during subsequent passes confirmed 
right side up capture. 

A simulation study was then made to determine tte 
performance of the spacecraft with the 50 ft boom 
in its actual orbit and to determine the boom 
length which would best meet the steady-state li- 
bration requirements . A run depicting the tran- 
sient motion of the spacecraft is given in Fig. 7- 
The steady -state response of the spacecraft, as 
mentioned earlier, is a function of the angle be- 
tween the line of nodes and the earth-sun line 
and the position of perigee - The node has a pre- 
cessional rate of about 2.24 degrees per day. 

Every 111 days the sun will have completed one 
revolution with respect to the orbit and there will 
have been four periods of increased roll libration 
amplitude . However, because of the increased alti- 
tude of apogee the effect of the thermal bending 
in roll is smaller than observed for the nominal 
orbit. An extrapolation indicated that the worse 
case occurs near day 280 in 1966 when perigee is 
in shade . At this time the simulation indicated 
that the deviation from local vertical is less thai 
5 degrees about 65 percent of the time and reached 
a peak of 7-5 degrees. The solar effect is mini- 
mized when the sun is in the orbit plane; at this 
time the peak libration angle is about 5*5 degrees. 
During the time of maximum solar effect, a boom 
length of 40 feet would result in deviations from 
the local vertical being under 5 degrees about 80 
percent of the time and a peak amplitude of 6 de- 
grees . Further reduction of the boom length would 
have little effect on the deviation from local 
vertical since the pitch libration dominates for 
the 40 ft length case . 

Post -Launch Attitude Detection System Performance 

While performance of the magnetometer system was 
satisfactory, the performance of the solar detec- 
tion system was disappointing. Soon after launch 
data indicated that both the +X and -X linear 
detectors were not functioning. In addition, all 
pairs of calibration detectors exhibited a one- 
step response in place of the expected two-step 
response curve illustrated in Fig. 3- The output csf 
the -X cosine detector also exhibited an abnormal 
response as a function of solar angle of incidence 
and the opening of a ground connection in this 
module on day 359 precluded the possibility of ob- 
taining meaningful attitude data after this date . 

The abnormal response of the -X cosine detector 
is illustrated in Figs . 9 and 10 . In Fig . 9 a plot 
of the 1M output illustrates that for some passes 
the response was a discontinuous function of time. 
Abnormal response of this detector was also observ- 
ed whenever the transition of illumination from 
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the +X to the -X cosine detector occurred. This 
phenomenon, illustrated in Fig. 10 for a pass on 
hay 329 when the sun is nearly in the plane of the 
orbit, produced abnormal attitude results . As il- 
lustrated by Figs . 10 and 11, the response of the 
Y and Z solar cosine detectors and all magnetome- 
ters indicated a continuous.motion of the satel- 
lite . The rate of change of the output of the +X 
and -X cosine detectors in the vicinity of the 
transition region compares favorably to the orbit- 
al rate of the satellite. A lag in the response in 
the -X cosine detector induced an apparent phase 
shift in the experimentally determined attitude of 
the satellite as illustrated in Fig. 12. For the 
portion of this pass over AEL and the beginning of 
the pass over Santiago, where the +X cosine detect- 
or is lit, good agreement is obtained between the 
experimental results and the digital simulation 
and linear eccentricity theory. For the remainder 
of the pass over Santiago the -X cosine detector is 
lit and the phase shift induced by this detector 
causes the difference between the experimental 
results and the simulation results to increase to 
approximately 7 to 8 degrees in hoth pitch and roll 
and to nearly 5 degrees in deviation from the lo- 
cal vertical. Abnormally high rates of motion are 
also associated with the changes in angles within 
this region. Statistical comparison of hoth por- 
tions of the pass with the digital simulation is 
given in Table II . Further indication of abnormal 
-X cosine detector performance was obtained from 
the magnitude of the post- launch calibration con- 
stants . Calibration was effected from data taken 
in the vicinity of day 340 from four stations: 
Santiago, Johannesburg, Canberra, and AH. A com- 
parison of the pre launch with the post- launch 
calibration constants is given in Table IIXA large 
decrease in magnitude of the calibration constants 
was observed for all cosine detectors with the 
exception of the -X detector. A decrease in the 
solar calibration constants has been observed in 
all satellites built by AEL, In GECS-A the -X de- 
tector calibration constant actually increased. 

All of the above factors support the conclusion 
that the -X cosine detector was not functioning 
properly and that the attitude data based upon 
this detector should he disregarded. 

Ccmparison of Theoretical and Experimental 
Attitude Data 

The agreement between the SAD experimental atti- 
tude and theoretical attitude of the spacecraft 
indicated by the simulation can be studied by 
examining Fig. 12, the combined AEL-Santiago pass 
on day 329. As discussed above, good agreement is 
obtained only over the portion of the pass during 
which the +X solar detector is lit. The SAD experi- 
mental data, which utilizes the malfunctioning -X 
solar detector, gives a physically absurd motion. 

An indication of the agreement between the ex- 
perimental and theoretical attitude can also be 
inferred from a comparison of the experimental and 
theoretical solar and geomagnetic data, see Figs. 

10 and 11. Small differences in the data over all 
of the double pass is observed indicating good 
three axis attitude agreement. This also confirms 


yaw capture of the spacecraft by the geomagnetic 
field. 

Since the Z detector data is sufficient to de- 
fine the orientation of the Z axis of the satel- 
lite, the very good agreement between the Z experi- 
mental data and the simulation results should be 
noted. At each data point the Z magnetometer and 
solar data independently define two cones of posi- 
tion for the Z axis of the satellite . The inter- 
section of these two cones yields, in general, 
two possible orientations of the sate lli tes' Z 
axis. One solution must be discarded by its indi- 
cation of an abnormally large deviation from the 
local vertical. The results of this simplified 
technique are presented and compared to the simu- 
lation results in Fig. 13. The spread at each data 
point is the uncertainty introduced by a one per- 
cent error in the experimental magnetometer and 
solar data - yielding nine possible combinations 
of input data. In general, the results of this 
simplified technique indicate extreme sensitivity 
to small errors in input data. A comparison of the 
SAD determination of attitude using redundant 
data with the simplified technique for the AH 
pass indicates the necessity for redundant data. 

It implies an error in Z detector data in excess 
of one percent. This comparison also indicates 
that the use of the simplified technique to deter- 
mine attitude over the Santiago pass is question- 
able. However, good agreement was obtained at the 
end of the Santiago pass especially for the total 
deviation from local vertical. 

The small difference between the total deviation 
from local vertical as computed from the linear 
eccentricity theory and the simulation indicates 
that the effect of all other disturbing torques 
should be small. This is because the sun Is nearly 
in the plane of the orbit; a geometric configura- 
tion which minimizes the roll disturbance . The 
simulation indicates a maximum total deviation of 
the Z axis from the local vertical of 4.5 degrees 
near the end of the Santiago pass. 

Results for additional passes over AFL, when the 
+X solar detector is illuminated, are given in 
Figs. 14 to 16 together with results of the digi- 
tal simulation and the linear eccentricity theory. 
Similar discrepancies at the beginning of each 
pass are noted for pitch and roll, hut there is 
reasonable agreement in the total deviation from 
local vertical when this occurs . A comparison of 
the experimental data over the four passes shows 
similarities from pass to pass indicating that a 
steady-state motion has been obtained and that 
the attitude determination system is stable. A 
statistical comparison between the SAD experimen- 
tal attitude data and the simulation results for 
these passes is given in Table IV. 

The agreement between simulation and experimen- 
tal attitude results indicates the successful per- 
formance of the GEOS -A stabilization system. The 
credibility of the digital simulation of the 
attitude motion of the gravity-gradient satellite 
has been supported. 
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TELEMENTRY VOLTAGE (volt*) 


50,000 



1. BOOM END MASS 
2 DOPPLER TRANSMITTERS 

3. DUAL OSCILLATOR 

4. SOLAR ATTITUDE DETECT- 
OR 

5 OPTICAL BEACON FLASH 
ASSY 

6 BROADBAND SPIRAL AN- 
TENNA 

7 RANGE AND RANGE RATE 
TRANSPONDER 


8. RANGE AND RANGE RATE 
ANTENNA 

9 LASSER REFLECTOR PANELS 
10. MEMORY 

U. FLIGHT PLUG ACCESS PAN- 
EL 

12. ANTENNA CONNECTOR AC- 
CESS PANEL 

13. SOLAR CELL ARRAY 

14. RANGE TRANSPONDER 


Fig. 1 GEOS-A SPACECRAFT 



Fig. 2 OUTPUT OF 4 GEOS-A LAMPS FLASHED SIMULTANEOUSLY 

BASED ON PROTOTYPE LAMP-REFLECTOR MEASUREMENTS. 
ANNA BEAM PATTERN IS SHOWN FOR COMPARISON. 



Fig. 3 GEOS-A SOLAR ATTITUDE DETECTOR RESPONSE 

curves, mi 


AMPLITUDE OF ROLL FREQUENCY COMPONENT OF 
FOURIER EXPANSION OF SOLAR FORCING FUNCTION 

=-^<fo »' n 1 G ** n W 



Fig. 4 AMPLITUDE OF ROLL ANGLE VERSUS ANGLE OF NODE 
FROM EARTH-SUN LINE GEOS-A IN NOMINAL ORBIT AND 
50 FOOT BOOM LENGTH 
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Fig. 5 AMPLITUDE OF ROLL ANGLE VERSUS ANGLE OF NODE 
FROM EARTH-SUN LINE, GEOS-A NOMINAL ORBIT 


THE NOMINAL INITIAL CONDITIONS ARE COMPUTED FROM 
LINEAR ECCENTRICITY THEORY AT 1642 UT ON DAY 317. 
NUMBERS REFER TO CHANCE IN RATE BY INDICATING 
EQUIVALENT LIBRATION AMPLITUDE IN DEGREES 

♦ INDICATES THAT MANEUVER PERFORMED AT THE NEXT 
HIGHER OR LOWER MINUTE RESULTED IN TUMBLING 



Fig. 6 GEOS-A INVERSION MANEUVER STUDY RESULTS. 

MAXIMUM AMPLITUDE AFTER INVERSION AS FUNCTION 
OF COAST TIME AND INITIAL ANGULAR RATE. 
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Fig. 7 GEOS-A TRANSIENT SIMULATION RUN 
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MAGNETIC FIELD (oersteds) 



TIME (second*) 


Fig. 9 GEOS-A TELEMETRY DATA ILLUSTRATING 

DISCONTINUOUS RESPONSE IN -X COSINE DETECTOR 



TIME (seconds) 


Fig. 10 GEOS-A SOLAR COSINE DETECTOR DATA FOR APL 
SANTIAGO PASS ON DAY 329 



| LINEAR ECCENTRICITY THEORY 

• SAD ATTITUDE DATA 
SIMULATION RESULTS 



+ X DATA 

- X DATA \ • 

.X- V 

. .. • , . . 



57600 56000 58400 58800 59200 

TIME (seconds) 


Fig. 11 GEOS-A MAGNETOMETER DATA FOR APL /SANTIAGO 
PASS ON DAY 329 


Fig. 12 GEOS-A ATTITUDE DATA AND SIMULATION RESULTS 
FOR APL/SANTIAGO PASS DAY 329 
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TOTAL 

(degrees) ROLL (degrees) PITCH (degrees ) TOTAL (degrees) ^ ^ ROLL (degree*) ( PITCH (degre**) 



57600 58000 58400 58800 59200 


TIME (seconds) 

Fig. 13 SIMULATION AND Z DETECTOR RESULTS FOR 
APL SANTIAGO PASS DAY 329 



51800 51900 52000 52100 52200 52300 52400 52500 52600 

TIME (seconds) 

Fig. M GEOS-A ATTITUDE DATA AND SIMULATION 
RESULTS FOR APL PASS ON DAY 338 



44800 44900 45000 45100 45200 45300 45400 45500 45600 45700 
TIME (seconds) 

Fig. 15 GEOS-A ATTITUDE DATA AND SIMULATION RESULTS FOR 
APL PASS ON DAY 339 



45100 45200 45300 45400 45500 45600 45700 45800 45900 


TIME (seconds) 

Fig. 16 GEOS-A ATTITUDE DATA AND SIMULATION RESULTS 
FOR APL PASS ON DAY 340 
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Table III 



Fi,. 17 IN-PLANE CONFIGURATION OF SPACECRAFT IN ORBIT 


Table I 

Satellite and Orbit Characteristics 

Satellite characteristics in launch configuration 

I x = 15.6 slug ft 2 

I = 17. It- slug ft? 

I z = 20.8 slug ft 2 

M = 386.6 lbs, mass 

Satellite characteristics for 5 0 ft configuration 

I x = 615.3 slug ft 2 

Iy = 617.0 SlUg ft 2 

I z = 20.8 slug ft 2 

A = 12,lt9- cm 2 , radiation area 
r = 175-3 cm, radiation moment arm 
c = 70,000 dyne-cm-sec, damping coefficient 
6 - 28.9 cm, thermal bending tip deflection 

P x ,Py,P z = -56.3, 159-, 250. pole-cm, dipole comp. 


Oybit characteristics 

actual 

nominal 

semi-major axis, km 

8073.5603 

7674.5410 

eccentricity 

.071922542 

.024 

inclination, deg 

59.371434 

59- 

ascn. node, deg 

40.256994 

variable 

arg. perigee, deg 

166.75751 

variable 

mean anomaly, deg 

0 

0 

rate of node, deg/da y 

-2.2457252 

-2.62 

rate of perigee, deg/day 

. 65242004 

.84 

time of perigee, day- year 

337-1965 

variable 

sec 

1812.4566 


Table 
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Statistical Comparison Between Attitude Data and 
Simulation Results for APL/Santiago Pass, Day 329 

RMS in degrees 
+ X data - X data 

Pitch 0.88 5.84 

Roll 1-55 4.74 


Analog Cosine Detector Calibration Constants 
Description 

of set +X -X +Y - Y - Z 

Prelaunch .437 .360 .492 .508 .444 

Post-launch .312 .398 - .375 -349 


Table IV 

Statistical Comparison Between Attitude Data and 
Simulation Results for APL Passes on Days 338,339 
and 340 

RMS in degrees 

Day 338 Day 339 Day 340 
Pitch 1.30 1.49 1.25 

Roll 1.76 1-79 1.98 

Total 0.80 1.20 .71 


Appendix A 

Steady-State Libration Induced by Orbit Eccentric- 
ity 

In this appendix the pitch (or in-plane) motion 
of a gravity-gradient satellite induced by orbit 
eccentricity is considered. For small values of 
eccentricities and small librational amplitudes, 
the pitch equation becomes decoupled from the roll 
and yaw equations. An approximate analytical solu- 
tion then becomes possible. 

Euler* s equation for the pitch motion of a satel- 
lite, depicted in Fig. 17, is given by 

I p (8 + $) = T g = ~ 3 (\ ~ I r ) 2 <t) (A.l) 

where T g represents the gravity- gradient torque. 

For an orbit of modest eccentricity 


9 = M + 2 e sin M + 0( e ) 
r/a = 1 - e cos M + 0(e 2 ) 


0 = - 2 g n 2 sin M + o(e 2 ) (A. 4) 

(a/r)^ = 1 + 3 e cos M + 0(c 2 ) (A. 5) 

Substituting Eqs. (A.1+) and (A. 5) into Eq. (A.l) 
yields 

<j>” + 3 O lj)= 2 e sin M + 0(<j) e ) + 0( e 2 ) (A. 6) 


M = nt j n = k/a^ 

a Ml r - I y )/I p 


and it has been assumed that 


0(|()) = 0(e) 

The steady-state solution to Eq. (A. 6) is 
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CT 


(B.7a) 

(B.7b) 


§ = “J sin M + 0(e) + 0(<|),e) (A.8) 

which yields the induced pitch motion as a function 
of mean anomaly. 

Appendix B 

Simplified Study of the Magnetically Anchored 
Damper 

In this appendix the damping and perturbing 
effects of the eddy current damper are examined. 

An analysis applicable to polar and near-polar 
satellites yields the damping rate and disturbing 
effect in terms of the damping coefficient and the 
inertia of the space craft. 

Extension of the gravity-gradient boom usually 
occurs only after the spacecraft has despun to 
near orbital rates. ®ie large inertia change 
which accompanies boom extension reduces this rate 
of rotation in inertial space to near zero about 
all axes except the longitudinal axis. Since the 
satellite in its stabilized position has an angular 
velocity in-plane of one cycle per orbital revolu- 
tion, pitch libration about this equilibrium posi- 
tion will result. For this reason a study of the 
pitch motion best characterizes the damping and 
perturbing effects of the damper 

The magnetic field for the satellite depicted in 
Fig. 17 which is assumed to move in a magnetically 
polar circular orbit, is given by 

H = H q (cos 0e 0 - 2 sin 0e r ) (B.l) 

The angle between the local vertical and the mag- 
netic field vector is given by 

y = - arctan cot 0) (B.2) 

If the magnet of the magnet assembly of the damper 
is collinear with the local field, the motion of 
the assembly relative to the conducting shell is 
given by 

Y - $ — o <i> ( B. 3 ) 

1+3 sin 0 

Experiments show that the torque exerted by the 
magnet assembly on the conducting shell is propor- 
tional to their relative velocity, c(y -(f); 
where c is the damping coefficient. 

For pitch motion of the satellite on a circular 
orbit, Euler r s equations of motion become 

I <b = T + T (B.4) 

p T g c 

where 

T g = | (I y - I r ) sin 2$ (B.5a) 

T c = c(y - <j>) (B.5b) 

For <[) small, Eq. (b. 4) becomes 


where 

a = (I - I )/l 
v r y" p 

( )• , M = nt 

The time constant, that is the time it takes 
for any disturbance to decrease to e"^ of its 
initial value, is given by 
nl 

t c = — ^ (orbital revolutions) (B.8) 

The time constant decreases and the disturbing 
torque increases as the ratio c/l increases. 

This means that any improvement ifi transient 
performance by increasing c/l degrades the steady- 
state performance. These results are valid for 
polar and near polar orbits. Digital simulation 
results indicate that Eq. (B. 8 ) is conservative. 
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